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Abstract 

We propose to use the different kinds of vacua of the gauge theories 
coupled to matter as a laboratory to test confinement ideas of pure Yang- 
Mills theories. In particular, the very poor overlap of the Wilson loop with 
the broken string states supports the 't Hooft and Mandelstam confine- 
ment criteria. However in the Z2 gauge-Higgs model we use as a guide we 
find that the condensation of monopoles and center vortices is a necessary, 
but not sufficient condition for confinement. 

1 Introduction 

Center vortices |l| and magnetic monopoles are widely believed to be the 
collective degrees of freedom responsible for the non-perturbative features of 
SU(N) Yang Mills theories, in particular for confinement. 

Condensation of magnetic monopoles implies a dual Meissner effect: the 
chromolclcctric field is expelled from the vacuum and gives the well-known phys- 
ical picture of confinement in terms of dual Abrikosov vortices which describe 
the flux tubes joining static sources. 

Center vortices are string-like excitations formed out of the center of the 
gauge group. They produce a very efficient disordering mechanism of the gauge 
configurations which could lead to area law decay of large Wilson loops. 

Although the above descriptions for almost all known models remains at 
a conjectural stageQ in the last few years many numerical lattice studies have 
given strong support for the relevant role played by center vortices and monopole 
condensation in confinement. For instance in the SU(2) gauge model in 3+1 

1 The only explicit, analytic, example of these mechanisms can be found in the Polyakov 
|3| proof of confinement of compact U(l) gauge model in 2+1 dimensions. 
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dimensions it has been observed both magnetic monopole condensation ||, [|, [| 
Q and the phenomenon of center dominance Q , namely the fact that the string 
tension obtained from center projected configurations in maximal center gauge 
agrees with the same quantity calculated in the full theory. 

Center dominance is verified in a trivial way in abclian theories, where the 
whole dynamics is described by center vortices. Here the distinction between 
confined and unconfined phase is related to the maximal size of the clusters 
of vortices: it can be shown that the confinement requires the presence of an 
infinite cluster . This agrees with an earlier observation of percolating center 
vortices in the cold phase of SU(2) gauge model at finite temperature For 
the sake of brevity we call vortex condensation the appearance of such a kind 
of infinite cluster. 

There are many open, intertwined, questions about the actual validity of 
the 't Hooft criteria of confinement, namely that the monopole and/or vortex 
condensations imply the decay with an area law of large Wilson loops. In par- 
ticular, is it possible to derive monopole condensation from vortex condensate 
or vice versa ? If not, are both condensations necessary for confinement? Arc 
they also sufficient? 

In this work we propose to answer these questions by studying gauge systems 
coupled to matter. Indeed it is worth noting that while center vortices and 
monopoles are generally studied in pure Yang Mills models, they are well defined 
also in gauge theories coupled to matter. In these models there are different 
vacua, distinguished by different entities which condense. If, for instance, the 
region where the vortices condense does not coincide with that with monopole 
condensation one can infer that these two properties are logically independent. 
This is precisely what happens in a 2+1 Z2 gauge-Higgs model || . A simple 
argument shows |L2| that the only vacua compatible with an area-law decay of 
the large Wilson loops are those in which both magnetic monopoles and center 
vortices condense. In other terms, the 't Hooft criteria are both necessary for 
confinement. 

Assuming that they are also sufficient can explain a surprising phenomenon 
observed in almost all coupled gauge systems studied up to now: although 
the potential between static sources flattens at large distances because of the 
screening produced by pair creation, this flattening (called string breaking) is 
invisible in the Wilson loop: it continues to obey an area law in full QCD jl3| 
even at distances where the static charges are completely screened. The point 
is that in gauge theories coupled to matter the basis of the operators has to be 
enlarged E^] in order to get a reliable estimate of the potential. In this way 
it has been observed the breaking of the confining string in Higgs models Jj"5[ 
and in QCD |l6). So the fact that large Wilson loops obey an area law even 
in coupled systems, as first suggested in Jf?]], may be considered as a further 
support to the usual plausibility arguments for the confinement mechanisms. 
Of course numerical experiments cannot give a true proof of the sufficiency of 
these criteria: it could well happen that Wilson loops of much larger size exhibit 
string breaking explicitly. 

In this work we study the 2+1 Z2 gauge-Higgs model to investigate a sub- 
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tier issue:by probing the different vacua characterised by infinite vortex and 
monopole clusters we study the effect of other condensates on confinement. In 
this model there are two vacua having center vortex and magnetic monopole 
condensates; one of them has also an electric condensate, i.e. the Higgs field 
has a vacuum expectation value different from zero. 

The former fulfils all the requirements of the confinement criteria, so an area 
law is expected; here we find that the Wilson loop obeys a perfect area law even 
at distances larger than five times the string breaking scale. 

The latter can be identified with the torn phase predicted in Ref.[[l7|: the 
Wilson loop decays with an area law below a given scale which varies very rapidly 
as a function of the couplings of the model and is unrelated to the string breaking 
scale. Above this threshold we observe a perimeter law: the infrared properties 
of this vacuum are indistinguishable from that of the perturbative unconfincd 
vacuum^. 

2 The model 

The action of a 3D Z2 gauge theory coupled to matter in a cubic lattice A can 
be written as 

S(J3 G , Pi) = -0i ViUiM ~PgY, U ^- W 

(ij) plaq. 

where both the link variable Uij = Ui and the matter field tpi take values ±1 

and U a = ILgD U ^ 

This model is self-dual: the Kramers- Wannier transformation maps the 
model into itself. Its partition function 

Z((3 G ,(3 I )= £ e -W°M (2) 
{<Pi=±l, u e = ±i} 

fulfils the functional equation 

Z(/3 G ,0i) = (sinh2/3 G smh2/3 / )l w Z(/3 7 ,/fc) (3) 

with j3 = — I log(tanh/3). 

The phase diagram of this model (see Fig. |l|) has been studied long ago j22j. 
There is an unconfined region surrounded by lines of phase transitions toward 
the Higgs phase and its dual. These lines are second order until they are near 
each other and the self-dual line, where first order transition occurs 

2 A similar phase has been reported in 4D S{7(2)-Higgs model Jis| 

3 For a more detailed description of the phase diagram see the poster presented by A.Rago 
at Lattice 2002. 
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Figure 1: Phase diagram for the 2+1 Z 2 gauge- Higgs system. The dashed 
(solid) lines denote the first-order (second-order) transitions. The dotted line is 
the self-dual line. 

2.1 Z 2 vortices and monopoles 

In this model the construction of center vortex configurations is straightforward: 
to each frustrated plaquette (i.e. U\j = —1) assign a vortex in the dual link. 
Since the product of the plaquettes belonging to any elementary cube is 1, center 
vortices form closed subgraphs of even coordination number. Thus a connected 
vortex subgraph contributes to a given Wilson loop W(C) only if an odd number 
of lines are linked to it: the Wilson loop is a vortex counter modulo 2. 

Magnetic monopoles can be defined exactly as in pure gauge model. They 
live on the dual lattice A. To create a monopolc in a site x € A, corresponding 
to the center of an elementary cube of A, it is sufficient to draw an arbitrary, 
continuous line 7(2;, y) joining x to to another monopole located in y (or to 
00) and flipping the sign of the coupling of the plaquettes crossed by 7. This 
flipping is generated by the non-local operator 

* 7 (i,y) = exp(-2/? G ^L/ n ) . (4) 

□ £7 

As a consequence, the flux across any closed surface with x inside and y outside 
is equal to -1: this monopole field \& 7 (x,2/) creates one unity of Z 2 flux joining 
x to y. Monopole condensation occurs when 

Urn (%(x,y))^0 . (5) 

\x-y\^oo 

Why this condensation should imply confinement? A useful piece of information 
comes from the Kramers- Wannier duality. One can easily show that under this 
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duality map the monopole correlator transforms as follows 

Thus the monopole condensation is associated to the dual Higgs phase, where 
the Z2 symmetry is spontaneously broken. 

The same transformation maps the Wilson loop W(C) associated to any 
closed curve C £ A into the corresponding 't Hooft loop W(C) of the dual 
phase: 

(W(C)) M = {W{C)) hr0a (7) 

with 

W(C) = exp(-2/3 G J2 ViUiiVo) > <9E = C (8) 

where E is an arbitrary surface bounded by C. W(C) creates an elementary Z2 
flux along C which manifests itself as a topological defect: the action of W(C) 
on an arbitrary configuration maps the product r\^, = Il^gc along any loop 
C £ A having an odd linking number with C into — r?g. 

2.2 A microscopic picture of confinement 

Confinement at 0g,Pi requires area law decay of (^(C))-^^. To understand 
this property at a microscopic level it is convenient to resort to the Fortuin 
Kasteleyn (FK) random cluster representation of the model. Starting from 
the obvious identity 

e 0I ^ u ^ =e f3 '(l-p + P 8 h(PxUlsy<Pv ) , p=l- e - 2 ^ (9) 

It is easy to perform explicitly the sum on the matter fields ip x , which yields 

Z{f3 G ,(3 I )=Y,^ npGUn Y,™u(G)v bG 2 c ° , v = j^—, (10) 

Ue GCA P 

The summation is over all spanning subgraphs GCA. bo is the number of links 
of G, called active bonds (which are the matter degrees of freedom replacing 
the ip's in this representation), and cq is the number of connected components, 
called FK clusters. w\j{G) is a projector on the subgraphs which are compatible 
with a given gauge configuration U = {Ue}: only those subgraphs are allowed 
for which no closed path within each FK cluster is linked to an elementary Z2 
flux |2^]. Put differently no frustration is permitted in G. In a sense, this 
is a microscopic realization of a sort of dual Meissner effect: the FK clusters 
behave like pieces of dual superconducting matter of type I and no Z 2 flux 
can go through the circuits of G. This is the only constraint generated by 
the interaction between matter and gauge fields. In the limit Pg ~ * 00 all 
the plaquettes have Uq = 1 (trivial gauge vacuum) hence the sum over the 
subgraphs is unconstrained and one gets the standard Ising model. Introducing 
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Figure 2: The different vacua of the Z2 gauge-Higgs model, characterised by 
different kinds of infinite FK or CV clusters (see Tab.|l|). 



a further projector z&c(G) on the space of configurations {G C A}, which takes 
the value whenever there is a FK cluster (at least) linked to C, and the value 
1 in all the other cases, yields the very useful identity 



(W(C)) 



number of compatible config. 
total number of config. 



(11) 



Assume that there are only FK clusters of finite size ( this is the case of the 
region V in Fig. ||). If C is much larger than the mean size of the clusters, 
the configurations contributing to W(C), i.e. those with wc{G) = 1, are char- 
acterised by the fact that there is no cluster linked to C . The weight of this 
class of configurations, when compared with the total ensemble G C A is clearly 
suppressed by a factor e~ ap ^ c \ where p(C) is the length of C. Thus in this 
phase the Wilson loop obeys a perimeter law. 

Conversely, a large Wilson loop decaying with an area law implies by ne- 
cessity the presence of an infinite FK cluster. In order to see whether this 
condition is also sufficient for confinement in the dual phase, note that accord- 
ing to Eq.(pT|) an infinite FK cluster fk^ gives a non- vanishing contribution 
to W(C) only if there is at least a simply connected surface £ bounded by C 
which is not pierced by the loops of the cluster. If the links of fkoc arc weakly 
correlated, one is led to argue that the weight of the configurations compatible 
with W(C) is suppressed by a factor e -CTa ( s ), where a(E) is the area of the 
minimal surface with = C. This leads to the area-law decay of (W(C)). 

A crucial assumption in the above argument is the weak correlation of the 
links of fkoo which describe at the microscopic level the monopole condensation. 
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Table 1 : Distribution of the infinite clusters in the phase diagram 



phase 


magnetic 


electric 


center 


dual 




condensate 


condensate 


vortices 


center vortices 


I 


yes 


no 


yes 


no 


II 


yes 


yes 


yes 


no 


III 


yes 


yes 


no 


yes 


IV 


no 


yes 


no 


yes 


V 


no 


no 


no 


no 



We shall see that there is a phase in which, though such a condensation occurs, 
this assumption is no longer true. Correspondingly we observe a violation of 
the area-law. 

Wc can now apply the same line of reasoning to the center vortices (CV). 
Finite CV clusters can link with the 't Hooft loop only along its perimeter. 
Therefore they contribute only to the perimeter term. An infinite CV cluster 
is necessary for the area-law: confinement implies both infinite center vortex 
subgraph and magnetic monopole condensation. In pure gauge theory these 
two requirements coincide, while in the coupled system the vacuum structure is 
more intriguing. 

2.3 The vacua 

In the coupled theory we have two kinds of dynamical subgraphs: CV or FK 
clusters in the dual lattice describe the gauge field degrees of freedom; FK 
clusters in the direct lattice describe the charged Higgs matter. 

In order to recognise an infinite cluster in practical simulations one has to 
look at the cluster size s as a function of the lattice volume V. We say that 
there is an infinite cluster whenever s oc V for large enough V. Straightforward 
numerical experiments show that they are distributed in the phase diagram 
according to Fig.|| and Tab.|l|. 

There are four kinds of infinite clusters: i) FK cluster in the dual lattice 
(magnetic condensate), ii) FK cluster in the direct lattice (electric condensate), 
Hi) CV cluster in the dual lattice, iv) CV cluster in the direct lattice (dual CV). 

In the region V there is no infinite cluster of any type. This is the perturba- 
tive, weak coupling vacuum, where large Wilson loops and 't Hooft loops obey 
a perimeter-law. 

The region IV is characterised by an electric condensate: it is the normal 
Higgs phase, where the Wilson loops decay with a perimeter-law, while the 't 
Hooft loops obey an area-law. 

Region III has a magnetic condensate but no large center vortices, then there 
is no confinement, as confirmed by numerical tests. This region is dual to the 
region II, where we shall see that the Wilson loops follow a perimeter-law decay. 
Because of duality relation ([?]) we can infer that also W(C) decays in the same 
way. 
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Figure 3: Results of Monte Carlo simulations of 6.0 ■ 10 5 sweep of update on 
cubic lattice of size 40 3 : (a):Plot of the rectangular Wilson W(R,T) loop and 
the other two operators defined in Eq.s (|lj) and ( |l3| ) at T = 17 as a function of 
R. (b):The static potential extracted by the lowest eigenvalue of the correlation 
matrix obtained for 6 different values of T (T = 15, • • • , 20) 

The regions I and II have both CV and FK infinite clusters in the dual 
lattice. 

The former is a normal confining phase: the potential extracted from an 
enlarged basis shows the expected string breaking (see Fig.|]b), while the Wilson 
loop obeys a perfect area-law even at a large scale (see Fig||). This has been 
checked up to distances of the order of five times the string breaking scale 
(compare Fig.||b and Fig.||). 

The latter is a phase with the simultaneous presence of infinite clusters in 
the direct and the dual lattices. The infinite cluster of A is associated to the 
condensation of the Higgs field <p x , while the infinite CV cluster of A is typical 
of a confining phase. However there is no confinement in the IR limit: the no- 
frustration constraint induces strong correlations among Z2 flux lines, because 
only an even number of them can pass through the closet paths within the FK 
clusters. 

As a matter of fact, large Wilson loops obey a perimeter-law decay, even if 
at intermediate distances an area-law seems recovered. A first example of this 
behavior is reported in ||. The cross-over scale varies rapidly as a function of 
the coupling constant and is unrelated to the string breaking scale. 

3 Numerical results 



As we have anticipated, even in this model, like in the other gauge systems 
coupled to matter analysed up to now, the Wilson loop fails to exhibit string 
breaking in the normal, confining phase, due to its poor (or vanishing) projection 
onto the screened potential. One has to use additional operators with good 
projection on the ground state, as first observed in Rcf.H]. The static potential 
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Figure 4: Improved square Wilson loops evaluated using duality relation (Q) and 
the projection on the largest FK cluster. The simulation has been performed 
on a 40 3 cubic lattice using 3.0 x 10 5 Monte Carlo configurations. 



V(R) and its excitations can be extracted from measurements of the matrix 
correlator, represented pictorially as 



C(R,T) 



V 



t 

T 



\ 



/ 



where Cn(R,T) is the rectangular Wilson loop (W(R,T)), the [/-shaped 
operator is 

C 12 (R,T) - (<p(0)U(0,-jT)U(-jT,-jT + kR)U(kR-jT,kR)<p(kR)) , (12) 

where j and k are two orthogonal unit vectors, and U(x,y) is a shorthand 
notation for a straight line of Ut connecting the sites x and y. The correlator is 
symmetric, C 12 (R,T) = C 2 i(R,T), and 

C 22 (i?,T) = {<p(0)U(0jT)ip(jT)ip(kR)U(kRjT + kR)ip(jT + kR)) . (13) 

Denoting by Ao < Ai the two eigenvalues, the ground state potential is defined 
as 



V(R) 



lim -logAo 

1 — t-OO 1 



(14) 
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In practical simulations, the limit T — > oo is not realized. In our case T was 
typically less that twenty lattice spacings when the the signal was lost in noise. 

The first point chosen in the present investigation is given by the couplings 
Pg = 0.75245 and fa = 0.16683 . This point in the phase diagram of Fig.| 
is in the phase I, fairly close to the deconfinement transition, where the model 
has properties very similar to the confining phase of pure gauge theory p3| . 
The estimates of the above three operators in a lattice of size 40 3 as functions 
of R are reported in Fig.||a and the static potential obtained the measure of 6 
different T in Fig.||b. We observe the typical string breaking phenomenon with 
a Wilson loop which seems to follow the area-law decay even in the region where 
the string is broken. In order to see better this property, we used a powerful 
algorithm based on Eq.(|Tl|) and already used in pure gauge theory p4|| . We 
also modified each configuration in the Monte Carlo ensemble by eliminating 
all the FK clusters not belonging to the largest cluster. It has been shown that 
this transformation does not change the value of the string tension but greatly 
reduce the noise of the measurements. The results are reported in Fig.^ : the 
square Wilson loops obey a perfect area-law even at distances five times the 
string breaking scale at this point. 

In order to investigate the effect of the electric condensate on the Wilson 
loop, we considered two different points in the phase diagram, one in the region 
I, at (3g — 0.6867 and (3j = 0.20, which is a normal confining phase. The other 
is in the region II, at {$q = 0.66 and /3j = 0.27 where besides the CV (and the 
FK) infinite cluster in the dual lattice there is also a percolating FK cluster in 
the direct lattice, associated to the "electric" condensation. These two points 
are chosen in such a way that the total size of the CV vortices are the same 
(within the statistical errors) in the two cases. It turned out that also the size 
of the maximal CV cluster is approximately the same. We also verified that this 
quantity scales with the volume for large enough lattices. 

The dramatic effect produced by the electric condensate is demonstrated by 
the completely different behaviors of the square Wilson loop in the two cases, 
as is evident in Fig]|. In the region I we see again the characteristic area-law 
of the normal confining phase, while the data of region II are compatible, for 
the Wilson squares of larger size, with a decay with the perimeter. This is the 



behavior expected for the torn phase 1 17 



4 Conclusions 

In this study we proposed to use the gauge theories coupled to charged matter 
to test the confinement ideas of 't Hooft and Mandelstam. In particular, in the 
case of the Z2 gauge-Higgs model, we gave a detailed microscopic description 
of the center vortex and monopole condensates. The analysis of the different 
vacua of this theory led us to conclude that 

• confinement requires that both magnetic monopoles and center vortices 
condense. 
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Figure 5: Square Wilson loops in the torn phase (circles) and in the normal, 
confining phase (crosses) resulting from two simulations in a cubic lattice of 
size 40 3 with 2.0 x 10 5 Monte Carlo configurations. The couplings are chosen 
in such a way that the density of center vortices is the same in the two cases. 
The straight line is a fit to a perimeter-law decay, while the dotted line fits an 
area-law decay with a = 0.1541(2). 

• the condensation of the former does not imply necessarily the condensation 
of the latter: there are vacua with a percolating FK cluster (i. c. the mi- 
croscopic description of the magnetic condensate) which is not associated 
to an infinite CV cluster. 

• plausibility arguments suggest that these two kinds of infinite clusters are 
responsible of the area-law decay of Wilson loop only if the links of these 
clusters are weekly correlated. 

• there are vacua where the condensation of the gauge degrees of freedom 
(monopolcs or center vortices) is associated with a condensate of the mat- 
ter field (the Higgs field in our example). This yields strong correlations 
among vortices and monopolcs so that the plausibility arguments for con- 
finement we alluded above are no longer justified. In fact we observed a 
pcrimctcr-law decay of large Wilson loops. This could be identified with 
the torn phase described in |jl7) . 
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